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Abstract 

In this paper we outline equipotential and equidensity methods of construct- 
ing axisymmetric models for galaxies. The former method defines equipotentials, 
from which the corresponding densities of the galaxy models can be obtained using 
Poisson's equation; the latter defines the equidensity surfaces of the galaxy models 
directly. 
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1 Introduction 

Most of the earlier models for galaxies are spherical and purely empirical, in order to fit the 
surface brightness of galaxies observed. Although spherical models can be used to simulate 
the surface brightness of some observed galaxies, galaxies are mostly not spherical. Thus 
we construct axisymmetric models in order to advance our overall understanding of galaxy 
formation and evolution. There are several different ways to construct axisymmetric 
models for galaxies. One of the basic ideas is to extend existing spherical models to more 
general axisymmetric forms. The methods of this extension can usually be classified into 
two groups. One is the equipotential method—, which was first introduced by Kuzmin' 111 ^ 
and developed by Toomre^, Miyamoto & Nagai^, SatohED, Kutuzov & Ossipkov' Mii ^, 
Evans^S, Jiang^, Ossipkov & Binney^S, Jiang & Moss^ and Jiang et al.®. This device 
is to define equipotentials from which the corresponding densities of galaxy models can 
be derived using Poisson's equation. The significance of this device is that it allows 
flattening to be achieved while maintaining a simple form of the potential. In particular, 
comparing the circular velocity profile generated from the potential of the model with the 
rotation curve of any very flattened galaxy, we can easily determine whether the model 
may be applied to describe the galaxy. An alternative is the equidensity method, which 
defines the equidensity surfaces of the galaxy models directly; this group includes the 
flattened 7-models^. The advantage of this method is that the projected surface densities 
of the models can be easily derived from their equidensity surfaces, allowing determination 
of whether these surface densities can be used to fit the surface brightness of observed 
galaxies. 
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2 The Equipotential Method 

Although the construction of galaxy models by the equipotential method has been out- 
lined by Kutuzov & Ossipkov^, we introduce below our rather different ideas about the 
method. Most of these ideas have been given by Jiang®. The equipotential method of 
constructing axisymmetric models for galaxies is currently being developed further. Up 
to now, the method can be roughly divided into three groups. The first is to use the 
axisymmetric radius of \/ R 2 + (a + \z\) 2 in place of the spherical radius r = R 2 + z 2 
in the potentials of the spherical models, where a > 0. This was originally introduced by 
KuzmintHEfl. It is not directly relevant to elliptical galaxies, and it is used to model discs. 
The physical significance of this device is that at the point (R, —\z\) below Kuzmin's disc, 
the potential of Kuzmin's models^ 

is identical with that of a point mass located at distance a above the centre of the disc. 
This implies that V 2 $(i?, z) = for z ^ 0. As a tends to zero, Kuzmin's models becomes 
point mass models. The second is to use the axisymmetric radius \J R 2 + (z/q) 2 in place 
of the spherical radius r of the potentials of spherical models, where q is the axial ratio, 
as in Binney's logarithmic models^. Evans's power-law models^ are given by a simple 
power-law potential 

*H--(^ F3 ,P>0. (2.2) 

where m is the axisymmetric radius given by m 2 = R 2 + (z/q) 2 , ip c is the potential 
at the origin and R c is the scale radius. In fact, Evans's power-law models can be re- 
garded as extensions of Plummer's spherical models^, since ( 12.21) degenerates to the 
spherical Plummer models when q = p = 1. The third is to use the axisymmetric radius 



R 2 + (V z 2 + c 2 + d) 2 in place of the spherical radius r of the potentials of the spherical 
models. This device was first introduced by Miyamoto and NagaiE3 for the Plummer 
potential, when they constructed the following potential-density pair 

HR, z) = - ° M = =, (2.3) 

^R 2 + (Vz 2 + c 2 + d) 2 

and 



= fc^M\ dR 2 + (d + 3y^+^) (d + Vz^T^) 2 
P{ }Z) ~\Ak J [R2 + (d+ Vz 2 + c 2 ) 2 fl 2 (z 2 + c 2 f/ 2 ' 

When c = 0, ( 12.31) reduces to Kuzimin's potential mentioned above. When d — 0, ( 12.31) 
reduces to the spherical potential of Plummer. 

The first and third devices described above have attracted the attention of many 
outstanding investigators. Thus several ideas of constructing axisymmetric potentials are 
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given and they are the development of these devices. A new potential-density pair can 
be obtained if (12 .ip is differentiated with respect to a 2 . Toomre's model n was derived by 
ToomreE2 by taking the (n — l)st derivative of (|2.ip /q with respect to a 2 . Similarly, if 
( 12. 3p is differentiated n times with respect to c 2 , Satoh's models^ can be obtained. 

Following Kuzmin's models mentioned above, Kuzmin & Kutuzov^ constructed a 
more general axisymmetric potential as follows 

C M 

$(R, z) = (2.5) 

(R 2 + z 2 + a 2 + c 2 + 2Va 2 c 2 + c 2 R 2 + a 2 z 2 ) 1 / 2 

This is a Stackel model. Similarly, Evans et al.^ obtained flattened isochrone models 
with the potential 

^/r, x X + aY + c 2 y 

^• Z> = - Y(X + aY + a*y ^ 

where X = ^a 2 c 2 + c 2 R 2 + a 2 z 2 and Y = Va 2 + c 2 + R 2 + z 2 + 2X. In (T23D and (M a 
and c are non-negative length scales. 

Recently, oblate Jaffe models for galaxies were given by Jiang 7 , using Miyamoto 
& Nagai's device. Then Jiang & Moss^ constructed a class of prolate Jaffe models, 
obtained by replacing the spherical radius r of the potentials of Jaffe's spherical models 



with a different axisymmetric radius of y (VR 2 + a 2 + b) 2 + z 2 , where a and b are two 
positive constants. This method is similar to that of Miyamoto & Nagai. Furthermore, 
Jiang et slM gave more general flattened Jaffe models with potentials 



GM / J(VW+^ + b) 2 + (V^T^+d) 2 

$(£> 2 , z ) = i n [ V = | , (2.7) 

! (VR 2 + a 2 + b) 2 + (Vz 2 + c 2 + d) 2 + rj 



rj 



using a device similar to that of Miyamoto & Nagai, that is replacing the spherical radius 
r of the potentials of Jaffe's models^ with a more general axisymmetric radius m, where 



,n = V (V-R 2 + a 2 + b) 2 + (V z 2 + c 2 + d) 2 , and a, 6, c and d are positive constants. 

Since the spherical 7-model ^ | 20 | 2:i | 2 ^ include Jaffe's spherical model, we can also define 
a class of axisymmetric 7-models, replacing the spherical radius r of the potential of the 
spherical 7-model with the more general axisymmetric radius m given above. 



3 The Equidensity Method 

The equidensity method directly defines the ellipsoidal equidensity surfaces of the galaxy 
models. For example, the axisymmetric radius ^ R 2 + (z/q) 2 can be used in place of the 
spherical radius r in the densities of the spherical models, where q is the axial ratio, as in 
Binney's logarithmic models^. The family of densities of the flattened 7- model' 3 ' is 
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from the density of the spherical 7-model mentioned in Section [2J, where M and a are 
total mass and scale radius, m is the same as in (12.21) . q is the axial ratio and < 7 < 3. 

The advantage of this method is that the projected surface densities of the models are 
easily derived. Thus we can easily determine whether the surface densities can be used 
to fit the surface brightness of the galaxies observed. However, this method can lead to 
mathematical difficulties in finding analytical expressions for the potentials. 
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